SWAT 98/212 
CERN-TH/98-398 



The Proton-Spin Crisis: another ABJ anomaly? 



G.M. Shore 

Department of Physics, University of Wales Swansea 
Singleton Park, Swansea, SA2 8PP, U.K. 

and 

Theory Division, CERN 
CH 1211 Geneva 23, Switzerland 



Contents: 

1. Introduction 

2. The First Moment Sum Rule for g\ 

3. The Parton Model and the 'Proton Spin' 

4. The CPV Method and Topological Charge Screening 

5. Experiment 

6. Semi-Inclusive Polarised DIS 



1. Introduction 

For a decade, the 'proton spin' problem - the anomalous suppression observed in the 
flavour singlet component of the first moment of the polarised proton structure function 
gP(x;Q 2 ) - has puzzled and intrigued theorists and experimentalists alike. The consequent 
research effort has indeed been impressive: to date, the original EMC paper [1] alone has 
nearly one thousand citations. 

In this lecture, we first give a brief review of the 'proton spin' problem from the 
standard viewpoint of the parton model. We explain why the 1988 observation by the 
EMC of a violation of the Ellis- Jaffe sum rule [2] for g\ was initially mis-interpreted in 
terms of quark spins and how the problem is resolved in the context of the full QCD 
parton model [3]. 
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We then describe an alternative, complementary approach to the description of deep 
inelastic scattering (DIS), the 'CPV method[4-6], which allows the problem to be viewed in 
a new light. From this perspective, the Ellis- Jaffe sum rule is simply seen to be equivalent 
to the OZI approximation for the forward proton matrix element of the flavour singlet axial 
current. The 'proton spin' problem is therefore one of understanding the origin of the OZI 
violation observed in this channel[7]. As such, it is one more addition to the collection of 
'CTa(I) problems' in QCD - phenomena whose interpretation depends on the presence of 
the ABJ axial anomaly [8] and the implicit relation with gluon topology (see e.g. ref.[9]). 
As we shall show, the observed suppression in the first moment of g\ is due to topological 
charge screening by the QCD vacuum itself, and a quantitative resolution in terms of an 
anomalous suppression of the slope of the gluon topological susceptibility is proposed. 

An immediate consequence of this explanation is that the suppression in g\ is in fact 
a target independent phenomenon, which would in principle be true for polarised DIS on 
any hadronic target. Not only is the 'proton spin' problem nothing to do with spin - it is 
not even a special property of the proton! To test this idea, we have proposed [10, 11] a set 
of semi-inclusive polarised DIS experiments, which could be performed at e.g. polarised 
HERA, and which would provide independent confirmation of the mechanism of topological 
charge screening by the QCD vacuum. 



2. The First Moment Sum Rule for g\ 

The structure function g\ is measured in polarised DIS experiments through the in- 
clusive processes up -> fj,X (EMC, SMC) or ep -> eX (SLAC, HERMES). See Fig. 1. 




Fig.l Inclusive polarised DIS scattering. 

The polarisation asymmetry of the cross section is expressed as 

dAa Ypl67r 2 a 2 v , ^(M 2 x 2 \ fn _ 

where the omitted terms include the second polarised structure function g\. The notation 

is standard: Q 2 = —q 2 and x = 2 ^ - are the Bjorken variables, y = ^ and Y p = 

According to standard theory, g\ is determined by the proton matrix element of two 
electromagnetic currents carrying a large spacelike momentum. The sum rule for the first 
moment (w.r.t. Bjorken x) of g\ is derived using the twist 2, spin 1 terms in the OPE for 
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the currents: 

J P (q)J a (-q) Q2 ~ oo 2e^^ [cf S (a s )(4 + + hof {a s )A^ (2.2) 

where C^ 5 and Cf are the Wilson coefficients and A a (a = 0, 3, 8) are the renormalised 
SU (3) flavour axial currents. The sum rule is therefore: 

T?(Q 2 ) = jf'cfa ^(*,Q 2 ) = ^< 5 (a 3 + Ja 8 ) + ^Cfa°(Q 2 ) (2.3) 

where the axial charges a 3 , a 8 and a°((5 2 ) are defined as the form factors in the forward 
proton matrix elements of the axial current, i.e. 

1 1 

(p,s\Al\p,s) = SfI -a 3 (p,s\Al\p,s) = s M ^=a 8 (p, s) = ^a°(Q 2 ) (2.4) 

Here, p M and s M = £47^7514 are the momentum and polarisation vector of the proton 
respectively. 

It is important in what follows to be precise about the renormalisation group properties 
of the operators and matrix elements which occur. The flavour non-singlet axial currents 
are not renormalised, but because of the {7a (1) anomaly the flavour singlet current A ^ is 
not conserved and therefore can, and does, have a non-trivial renormalisation. Defining 
the bare operators by A°^ B = J2Ql^T a q and Qb = f^^trG^G^, (where T a ^° are 
SU(rif) flavour generators, T° = 1, and <i-symbols are defined by {T a ,T 6 } = d a b c T c ), the 
renormalised operators A^, Q, are given by[12] 

Af° = A%° Al = ZA% Q = Q B -^-(l-Z)d»A% (2.5) 

where Z is a divergent renormalisation constant, the associated anomalous dimension being 
denoted by 7. Q is the gluon topological charge density - for field configurations which 
tend to a pure gauge at infinity, it satisfies 



/ 



d 4 x Q = n e Z (2.6) 



where n is the topological winding number, or 'instanton' number. 

The link between quark dynamics in the flavour singlet {7a (1) channel and gluon 
topology is provided by the famous {7a(1) axial (or ABJ) anomaly^: 

d^Al - 2n f Q ~ (2.7) 



( 1 ) In this lecture, we will quote results only in the chiral limit of QCD, i.e. with massless quarks. The 
inclusion of quark masses, which produces only a small change in our final predictions, is described in 
detail in ref. [13]. 
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where the symbol ~ denotes weak operator equivalence. Notice that with the definitions 
(2.5), this condition is the same expressed in terms of either the bare or renormalised 
operators. 

A more complete formulation of the anomaly is given by introducing the generating 
functional W[S®,6, S a ] of connected Green functions (correlation functions) of the 
axial currents and the pseudoscalar operators Q and (/>?■, where 4>^ B = q.l$T a q (together 
with the corresponding scalar (p a ). Here, S®,9,S§,S a are the sources for A a ^ Q, <p^, <p a 
respectively. The (anomalous) Ward identities are then expressed as 

d.JS* - ^ f S a o W + d adc S«— - d adc Si— = (2.8) 

where the final two terms account for the chiral variations of the operators <f>%, (fi a . (For 
a more extensive description of this formalism, see refs.[5,9].) Ward identities for 2-point 
Green functions (composite operator propagators) are then derived by taking functional 
derivatives w.r.t. the relevant sources. For example, 

d^W s o e -2n f W ee = (2.9) 

where for simplicity we adopt a notation where functional derivatives are denoted by sub- 
scripts and integrals over spacetime are assumed as appropriate. Since there is no massless 
Ua(X) Goldstone boson in the physical spectrum (resolution of the Ua{X) problem), the 
first term vanishes at zero momentum, leaving^ 2 ) ( 3 ) 

W ee \ k=0 = (2.10) 

The renormalisation group equations (RGEs) implied by eq.(2.5) can also be conve- 
niently written in a functional form. The fundamental RGE for the generating functional 
W is[5,9] 

VW = l(sl- ^-d^Wso + ^(S%W S « + S a W Sa ) +... (2.11) 
where V = [a-§- + 0jr- ) and 7^ is the anomalous dimension corresponding to the 

V °^ ag J V,8,S 5 ,S 

pseudoscalar (or scalar) composite field renormalisation. The notation + . . . refers to the 
additional terms which are required to produce the contact terms in the RGEs for Green 
functions involving more than one composite operator. These vanish at zero momentum 
and need not be considered here. 

The RGEs for 2-point Green functions follow immediately by functional differentiation 
of eq.(2.11). For example, differentiating twice w.r.t. 6>, and using the Ward identity (2.9), 
we find 

VW ee = 2-fW ed + ... (2.12) 



( 2 ) Here, and elsewhere throughout the text, I have included a discussion of points omitted during the 
lecture itself which were subsequently raised in the discussion sessions. 

^ In fact, the same conclusion would hold even if there were a physical massless boson coupling deriva- 
tively to the J/a(1) current since, in momentum space, the residue of the pole at k 2 = in the first term 
in eq.(2.9) due to this boson would itself be 0(k 4 ). 
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We are focusing on results for Wgg here because it is a key correlation function in QCD 
which will play a central role in our analysis of the 'proton spin' problem. It is called the 
topological susceptibility, and is usually denoted by x- I n conventional notation, 

x (k 2 ) = Wee{k 2 ) =i j ' d 4 x e* kx (0\T Q(x) Q(0)]0) (2.13) 

Returning to the structure function sum rule (2.3), we can now deduce the RG be- 
haviour of the axial charges. From either (2.5) or (2.11), we see immediately that 

|a 3 ' 8 = |a°(Q 2 ) = 7 a°(Q 2 ) (2.14) 

where t = lnQ 2 /A 2 . The singlet axial charge is therefore scale dependent. This is crucial 
in understanding the 'proton spin' problem in QCD. 

The flavour non-singlet axial charges a 3 and a 8 are known in terms of the F and D 
constants found from neutron and hyperon beta decay: 

a 3 = F + D a 8 = 3F-D (2.15) 

The interest of the sum rule therefore centres on the flavour singlet axial charge a°(Q 2 ). 
In the absence of an alternative experimental determination of a (Q 2 ), the simplest ansatz 
is to assume that it obeys the OZI rule, i.e. a°(Q 2 ) = a 8 . In precise terms, the OZI limit 
of QCD is defined [14] as the truncation of full QCD in which non-planar and quark- loop 
diagrams are retained, but diagrams in which the external currents are attached to distinct 
quark loops, so that there are purely gluonic intermediate states, are omitted. (This last 
fact makes the connection with the familiar phenomenological form of the OZI, or Zweig, 
rule.) This is a more accurate approximation to full QCD than either the leading large 
1/N C limit, the quenched approximation (small n/ at fixed N c ) or the leading topological 
expansion (N c — > oo at fixed nf/N c ). In the OZI limit, the {7a (1) anomaly is absent, there 
is no meson-glueball mixing, and there is an extra {7a(1) Goldstone boson. 





(b) 

Fig. 2 'Connected' (a) and 'disconnected' (b) contributions to the matrix element Diagrams of 

type (b) are suppressed in the OZI limit. 

Applied to the axial charges, the OZI rule states that of the two valence quark diagrams 
shown in Fig. 2 describing the coupling of the axial current to the proton, the diagram in 
Fig. 2(b) is suppressed. Since this is the only way the s-quark component of either A ^ or 
A 8 ^ can couple to the proton (uud), the OZI prediction aP = a 8 follows immediately. 
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Inserted into the sum rule (2.3), the OZI approximation a°(Q 2 ) = a 8 gives a theoretical 
prediction for which is known as the Ellis-Jaffe sum rule. This is now known to be 
violated, with a°(<5 2 ) strongly suppressed (by a factor of around 0.5) relative to a 8 . This 
experimental result, which was first discovered by the EMC collaboration and confirmed 
and improved in the subsequent decade of experimental work, is what has come to be 
known as the 'proton spin' problem (or even crisis!). 

In fact, it is not at all surprising that the OZI rule should fail in this case[7,4,5]. 
The first clue is the anomaly-induced scale dependence of a°(Q 2 ). If the OZI rule were 
to hold, at what scale should it be applied? Moreover, it is known that the pseudovector 
and pseudoscalar channels are linked through the Goldberger-Treiman relations. Since 
large anomaly-induced OZI violations are known to be present in the pseudoscalar channel 
(£/a(1) problem, r{ mass, etc.) it is natural to find them also for a°(<5 2 ) in the pseudovector 
channel. It is also immediately clear from its scale-dependence that a°(Q 2 ) cannot really 
measure spin. 

While this immediately resolves the 'proton spin' problem, clearly we want to under- 
stand the origin of the suppression in a°(Q 2 ) much more deeply. In the next two sections, 
we describe two complementary approaches to this question - the conventional QCD parton 
model and the CPV method developed in refs.[4-6]. 



3. The Parton Model and the 'Proton Spin' 

In the simplest form of the parton model, the proton structure for large Q 2 DIS is 
described by parton distributions corresponding to free valence quarks only. The polarised 
structure function is given by 

n f 

where Aqi(x) = qj r (x)+q^ r (x)—q~ (x)—q~ (x) is defined as the difference of the distributions 
of quarks (and antiquarks) with helicities parallel and antiparallel to the nucleon spin. It 
is convenient to work with the flavour non-singlet and singlet combinations: 

Aq NS (x) = ~ l ) A ^ = E ( 3 - 2 ) 

i=l \ e ' i=i 

In this model, the first moment of the singlet quark distribution Aq s = J* dx Aq s (x) can 
be identified as the sum of the helicities of the quarks. Interpreting the structure function 
data in this model then leads to the conclusion that the quarks carry only a small fraction 
of the spin of the proton - the 'proton spin' problem. There is indeed a real contradiction 
between the experimental data and the free valence quark parton model. 

However, this simple model leaves out many important features of QCD, the most 
important being gluons, RG scale dependence and the chiral {7a (1) anomaly. When these 
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effects are included, in the QCD parton model, the naive identification of Aq s with spin 
no longer holds and the experimental results for g\ are readily accommodated. ( 4 ) 




Fig. 3 QCD parton model interpretation of DIS. The single lines are partons, which may be quarks or 
gluons. 

The QCD parton model picture of DIS is shown in Fig. 3. The polarised structure 
function is written in terms of both quark and gluon distributions [3] as follows: 
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y 



c ns Q AqNS{yj t) + c s (^)A q s (y, t) + C° {^)Ag(y, t)] (3.3) 



where C s \C 9 and C NS are perturbatively calculable functions related to the Wilson co- 
efficients in sect. 2 and the quark and gluon distributions have a priori a t = lnQ 2 /A 2 
dependence. 

The RG evolution (DGLAP) equations for these polarised distributions are: 



5 A ^ = S 



1 ^y_pNs 
y qq 



(3.4) 



and, 



±(Aq s {x,t)\_a.[ 1 dy(Pf g {Z) P q9 {% ) \( A<T%, t) \ ( . 

dt \ Ag(x,t) ) 2tt J x y {P gq (V P M (f) J { Ag(y,t) ) 



Of course there is a separate angular momentum sum rule for the proton. The details are given in 
e.g. ref. [15]. The proton spin is given by the matrix element of the angular momentum operator J: 



d 3 xM 0jk 



where is the energy-momentum tensor. This admits (up to equation of motion terms) a gauge-invariant 
decomposition into three terms which look like quark spin, quark orbital and total gluon angular momenta: 

We recognise here the operator A^ B = ^<?7^75<7 which, for free fields, can be identified as the quark 
helicity operator. However, the composite operators in this expression must be renormalised and their scale 
dependence identified, and the above identifications then become more problematic. While it is interesting 
to pursue both the theory and experimental consequences of this sum rule [15], it is obvious that the axial 
charge a°(Q 2 ), given by the matrix element of the renormalised current is not measuring the spin of 
the proton. 
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showing the mixing between the singlet quark and the gluon distributions. The splitting 
functions P are also calculable in perturbative QCD, their moments being related to the 
anomalous dimensions of the series of increasing spin operators appearing in the OPE 
(2.2). 

In this language, the first moment sum rule for g\ reads: 

T\{Q 2 ) = 1 [c? S Aq NS + C?Aq s + C(Ag\ (3.6) 

where Aq NS , Aq s and Ag are the first moments of the above distributions. Comparing 
with eq.(2.3), we see that the axial charge a°(<5 2 ) is identified with a linear combination 
of the first moments of the singlet quark and gluon distributions. It is often, though not 
always, the case that the moments of parton distributions can be identified in one-to-one 
correspondence with the matrix elements of local operators. The polarised first moments 
are special in that two parton distributions correspond to the same local operator. 

The RG equations for the first moments of the parton distributions follow immediately 
from eqs. (3.4,3.5) and depend on the matrix of anomalous dimensions for the lowest spin, 
twist 2 operators. This introduces a renormalisation scheme ambiguity. The issue of 
scheme dependence has been studied thoroughly by Ball, Forte and Ridolfi[16] and an 
excellent summary can be found in ref.[17]. It is shown there that it is possible to choose 
a scheme known as the Adler-Bardeen or AB scheme (strictly, a class of schemes [16]) for 
which the parton distributions satisfy the following RG equations: 

| V" = ±A<f = 

with the implication Cf = — n/^jCf . It is then possible to make the following identifica- 
tions with the axial charges: 

a 3 = Au - Ad 
a 8 = Au + Ad- 2As 
a°(Q 2 ) = Aq s -n f ^-Ag(Q 2 ) (3.8) 

where Au = L dx(Au(x, t) + A-u(x,t)) etc. Notice that in the AB scheme, the singlet 
quark distribution Aq s = Au+Ad+As (which is often written as AE) is scale independent. 
All the scale dependence of the axial charge a°(Q 2 ) is assigned to the gluon distribution 
Ag(Q 2 ). We emphasise that eq.(3.7) is true only in the AB renormalisation scheme and 
that it is only in this scheme that the identifications (3.8) hold. 

This was the identification originally introduced for the first moments by Altarelli 
and Ross [3], and resolves the 'proton spin' problem in the context of the QCD parton 
model. In this picture, the Ellis-Jaffe sum rule follows from the assumption that in the 
proton both As and Ag(Q 2 ) are zero. This is the natural assumption in the context of the 
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free valence quark model. It is equivalent to the naive OZI approximation a°(Q 2 ) = a 8 
described above. However, in the full QCD parton model, there is no reason why Ag(Q 2 ), 
or even As, need be zero in the proton. Moreover, given the RG scale dependence of 
a°(<5 2 ), this assumption is in contradiction with QCD where the anomaly requires a°(Q 2 ) 
to scale with the anomalous dimension 7. 

Since neither AE nor Ag(Q 2 ) are currently measurable in other processes, the parton 
model is unable to make a quantitative prediction for the first moment T^. While the 
model can accomodate the observed suppression, it cannot yet predict it. 

An interesting conjecture, proposed in the original paper of Altarelli and Ross [3], is 
that the observed suppression in a°(Q 2 ) is due overwhelmingly to the gluon distribution 
Ag(Q 2 ). If so, the strange quark distribution As ~ in the proton and so AE ~ a 8 . 
Although by no means a necessary consequence of QCD, this is entirely plausible because 
it is the anomaly (which is due to the gluons and is responsible for OZI violations) which 
is responsible for the scale dependence in a°(Q 2 ) and Ag(Q 2 ) whereas (in the AB scheme) 
AE is scale invariant. The essence of this conjecture will reappear in the next section 
where we describe the CPV method. 

To test this conjecture, we need to find a way to measure Ag(Q 2 ) itself, rather than 
the combination a°(Q 2 ). One possibility[16,18] is to exploit the different scaling behaviours 
of Aq s (x) and Ag(x,Q 2 ) to distinguish their contributions in measurements of gf(x,Q 2 ) 
at different values of Q 2 . A second is to extract Ag(x,Q 2 ) from processes such as open 
charm production, ^*g — > cc, which will be studied in various forthcoming experiments at 
COMPASS, RHIC, etc. 



4. The CPV Method and Topological Charge Screening 

In this section, we shall discuss a less conventional approach to DIS based on a 
decomposition of matrix elements into products of composite operator propagators and 
their associated 1PI vertex functionals. This formalism has been developed in a series of 
papers[4-6,13] on the 'proton spin' problem. The starting point, as indicated above, is the 
use of the OPE in the proton matrix element of two currents. This gives the standard 
form for a generic structure function moment: 

f dx x^F^Q 2 ) = J2c?(Q 2 )(p\0?mP) (4-1) 
Jo 



where are the set of lowest twist, spin n operators in the OPE and C^tQ 2 ) the corre- 
sponding Wilson coefficients. In the CPV approach, we now factorise the matrix element 
into the product of composite operator propagators and vertex functions, as illustrated in 
Fig. 4. 
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Fig. 4 CPV description of DIS. The double line denotes the composite operator propagator and the lower 
blob the 1PI vertex. 

To do this, we first select a set of composite operators 0\ appropriate to the physical 
situation and define vertices as 1PI with respect to this set. Formally, this is achieved 

by introducing sources for these operators in the QCD generating functional W, then 
performing a (partial) Legendre transform[5] to obtain a generating functional r[Cy . The 
1PI vertices are the functional derivatives of r[d);]. The generic structure function sum 
rule (4.1) then takes the form 

i 3 

= EEw ( 4 - 2 ) 

i 3 

in a symbolic notation. 

This decomposition splits the structure function into three pieces - first, the Wilson 
coefficients Cj(Q 2 ) which control the Q 2 dependence and can be calculated in pertur- 
bative QCD; second, non-perturbative but target-independent QCD correlation functions 
(0|TOj n) &i\0); and third, non-perturbative, target-dependent vertex functions Tq de- 
scribing the coupling of the target proton to the composite operators of interest. The vertex 
functions cannot be calculated directly from first principles. They encode the information 
on the nature of the proton state and play an analogous role to the parton distributions 
in the more conventional parton picture. 

It is important to recognise that this decomposition of the matrix elements into prod- 
ucts of propagators and proper vertices is exact, independent of the choice of the set of 
operators Oi. In particular, it is not necessary for Oi to be in any sense a complete set. 
All that happens if a different choice is made is that the vertices themselves change, 

becoming 1PI with respect to a different set of composite fields. Of course, while any set of 
Oi may be chosen, some will be more convenient than others. Clearly, the set of operators 
should be as small as possible while still capturing the essential physics (i.e. they should 
encompass the relevant degrees of freedom) and indeed a good choice can result in vertices 
^ 6 pp w hich are both RG invariant and closely related to low energy physical couplings, 
such as g^NN- In this case, eq.(4.2) provides a rigorous relation between high Q 2 DIS and 
low-energy meson-nucleon scattering. 

For the first moment sum rule for g\, it is most convenient to use the ?7a(1) anomaly 
equation immediately to re-express a°(Q 2 ) in terms of the forward matrix element of the 
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topological charge Q, i.e. 



a°(Q 2 ) 



1 



2M 



2n f (p\Q\p) 



(4.3) 



where M is the nucleon mass. 

Our set of operators Oi is then chosen to be the renormalised flavour singlet pseu- 
doscalars Q and $5, where $5 is simply the operator 0° of section 2 with a special, and 
crucial, normalisation. The normalisation factor is chosen such that in the OZI limit of 
QCD, where the anomaly is absent, $5 would have the correct normalisation to couple 
with unit decay constant to the {7a (1) Goldstone boson which would exist in this limit. 
This also ensures that the vertex is RG scale independent. (The proof may be found in 
refs.[5,13].) The vertices are defined from the generating functional T[S£, Q, 0g, (fi a ] where 

r[5-, Q, C <p a ] = w[s a ^ e, s? , s a ] - f dx (oq + + s a <p a ) (4.4) 



We then have 



1 1 



singlet g 



2n f Cf{a s ) 



(o\t q Q\o)f Qpp + (o\t q <s> 5 \o)f<, 



hVV 



(4.5) 



where the propagators are at zero momentum and the vertices are 1PI w.r.t. Q and $5 
only. For simplicity, we have also introduced the notation iuYq pp u = YQppU'-y^u, etc. This 
is illustrated in Fig. 5. 



XQ 



XQ 



Q 





Fig. 5 CPV decomposition of the matrix element (p\Q\p). The propagator in the first diagram is x(0); 
in the second it is y/ x'(0). 

The composite operator propagator in the first term is simply the (zero-momentum) 
QCD topological susceptibility x(0) which, as we have seen in section 2, vanishes for QCD 
with massless quarks. Furthermore, with the normalisation specified above for $5, the 
propagator (0\T Q $s|0) at zero momentum is simply the square root of the slope of the 
topological susceptibility. 

To see this, notice that by virtue of their definition in terms of the generating functional 
(4.4), the matrix of 2-point vertices in the pseudoscalar sector is simply the inverse of the 
corresponding matrix of pseudoscalar propagators, i.e. 



QQ 

0SQ 



This implies 



r, 



Wee W es c 



W s o e 



Wgo go 

5 5 



-Wee{det W ss )' 



(4.6) 



(4.7) 
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letting S represent the set {6>, S®}. Differentiating w.r.t. k 2 and taking the limit k 2 = 0, 
exploiting the fact that W<9#(0) vanishes, we find 

^^o = x'(0)W- s l(0) (4.8) 

Finally, normalising the field $5 proportional to (p® such that ^3T$ 5 $ 5 \ k _ Q = 1 (see 
refs.[5,13] for further discussion), we find the required relation 

(0|TQ$ 5 |0)| fc=0 = v^(0) (4.9) 

We therefore find: 

H! singlet = l^M 2n f t^ 5PP (4.10) 

The slope of the topological susceptibility x'(0) is not RG invariant but, as shown in 
eq.(2.12), scales with the anomalous dimension 27, i.e. 

j t V7W) = iV7W) (4.11) 

On the other hand, the proper vertex has been chosen specifically so as to be RG invari- 
ant[5,13]. The renormalisation group properties of this decomposition are crucial to our 
proposed resolution of the 'proton spin' problem. 

Our proposal[4,5] is that we should expect the source of OZI violations to lie in the 
RG non-invariant, and therefore anomaly-sensitive, terms, i.e. in x'(0) rather than in the 
RG invariant vertex. This is the key assumption that allows us to make a quantitative 
prediction for on the basis of a calculation of the topological susceptibility alone. The 
RG invariance of the vertices is a necessary condition for this assumption to be reasonable. 
Further phenomenological evidence from Ua(X) current algebra supporting this conjecture 
is discussed in refs.[5,19]. Notice that we are using RG non-invariance, i.e. dependence on 
the anomalous dimension 7, merely as an indicator of which quantities are sensitive to the 
anomaly and therefore likely to show OZI violations. Since the anomalous suppression in 
is thus assigned to the composite operator propagator rather than the proper vertex, 
the suppression is a target independent property of QCD related to the anomaly, not a 
special property of the proton structure. 

In this picture, therefore, the basic dynamical mechanism responsible for the sup- 
pression of the 'proton spin' can be identified as topological charge screening by the QCD 
vacuum. That is, when a matrix element of the topological charge is measured, the QCD 
vacuum screens the topological charge through the zero or anomalously small values of 
the susceptibility %(0) and its slope x'(0) respectively (see Fig. 4). The mechanism is 
analogous to the screening of electric charge in QED. There, because of the gauge Ward 
identity, the screening is given entirely by the ('target independent') dressing of the photon 
propagator by vacuum polarisation diagrams, leading to the relation = eB\/~Zz (with 
Z3 < 1) between the renormalised and bare charges, in direct analogy to eq.(4.12) below 
with the topological susceptibility playing the role of the photon propagator. 
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To convert this into a quantitative prediction we use the OZI approximation^ 5 ) for the 
vertex T$ 5Pp . In terms of a similarly normalised octet field $|, this is T® 5PP = y/2T^s pp . 
The normalisation is crucial in allowing the use of the OZI relation here. The corresponding 
OZI prediction for a/ x'(0) would be f^/VS- These OZI values are determined by compar- 
ing the result (4.10) (at least that part relating to the proton matrix element) with the 
conventional Goldberger-Treiman relation for the flavour octet axial charge in the chiral 
limit (see refs.[5,13]). This gives our key formula^ 

for the flavour singlet axial charge. Incorporating this into the formula for the first moment 
of the polarised structure function, we find 

r? sin 9 let = lc?M a 8 ^ VtfW) (4-13) 

The final step involves an explicit calculation of %'(0). This was done in ref.[6] using 
the QCD spectral sum rule (QSSR) method. We find[6,13] 

\/x'(0) = (26-4 ± 4.1) MeV (4.14) 

Substituting this in eq.(4.13), and also using the QSSR prediction for f n in order to 
minimise systematic errors, we arrive at our prediction for the first moment of the polarised 
structure function in the chiral limit: 

a°(Q 2 = lOGeV 2 ) = 0.33 ± 0.05 

T\(Q 2 = lOGeV 2 ) = 0.144 ± 0.009 (4.15) 



( 5 ) Without this approximation, our result (4.12) reads: 

°°W 2 ) ( Ve r— \ ( i*,„ \ _ 
— -{uW>)\v) {-^f—J - «-v 

A (target-independent) propagator suppression factor sp < 1 is a signal of topological charge screening; 
a (target-dependent) vertex suppression factor sy < 1 would be a sign that the proton coupling to the 
flavour singlet Goldstone boson Vozi m ^ e limit of QCD is small, equivalent in parton terms to 
a significant (negatively polarised) strange quark content in the proton. Of course both effects may be 
present. However, our conjecture is that the 'proton spin' problem is best explained by sp < 1, sy — 1- 

( 6 ) It is interesting to compare eq.(4.12) with the well-known Witten-Veneziano formula[20,21] which, at 
leading order in the l/N c expansion, relates the mass of the rf (in massless QCD) to x YM > tne topological 
susceptibility of pure gluodynamics: 

V6 



In 
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This can be immediately compared with the Ellis- Jaffe sum rule prediction 

a°| EJ = 0.58 ±0.03 (4.16) 

We therefore find a suppression of a°(Q 2 ) relative to the OZI expectation by a factor of 
around 0.55. 

A complete description of our calculation of the slope of the topological susceptibility 
can be found in refs.[6,13]. The QSSR analysis of x'(0) in the chiral limit of QCD is 
essentially straightforward and shows a clear range of stability with respect to the Laplace 
sum rule parameter r, indicating that the prediction is reliable. The stabilisation scale, 
r" 1 ~2.5 — 5 GeV 2 , is sufficiently big for higher dimensional condensates to be suppressed, 
and the calculation displays the hierarchy of gluonium to light meson hadronic scales 
anticipated by ref.[22].( 7 ) 

Lattice gauge theory methods may also be used to calculate the topological sus- 
ceptibility. However, x'(0) is a particularly difficult correlation function to calculate 
on the lattice, requiring algorithms that implement topologically non-trivial configura- 
tions in a sufficiently fast and efficient way. Very preliminary results from the Pisa 
group [26] of calculations in full QCD with dynamical quarks indicate a value of the order 
^/x'(0) ~ (19 ±4) MeV. Given the preliminary nature of the lattice simulations, the rough 
agreement with the QSSR result is encouraging. Qualitative explanations of topological 
charge screening and the anomalously small value of x'(0) in QCD may also be given using 
models of the instanton vacuum, e.g. the instanton liquid model of ref.[27]. 

Finally, to complete this discussion, it is useful to recognise the complementary nature 
of the QCD parton model and the CPV method presented here. Both involve at present 
incalculable non-perturbative functions describing the proton state - the quark and gluon 
distributions in the parton picture and the 1PI vertices in the CPV method. Both exhibit 
a degree of universality - the same parton distributions may be used in different QCD 
processes such as DIS or hadron-hadron collisions, while the vertices (when they can be 
identified with low-energy couplings) also provide a link between high Q 2 DIS and soft 
meson- nucleon interactions. 

The principal attribute of the parton model is that it allows a detailed description of 
the structure of the proton in terms of its quark and gluon constituents. On the other hand, 

Nevertheless, our application of QSSR to the Ua(1) sector of QCD has been criticised repeatedly 
by Ioffe (see e.g. refs. [23,24]) on two grounds: (i) that there are important neglected contributions from 
'instantons', i.e. higher dimensional condensates, and (ii) that when the strange quark mass is included, 
QSSR results for current-current correlators show quite unrealistic SU(3) breaking. Neither criticism 
is valid, and both are refuted in detail in our recent paper[13]. We have commented in the text how 
the size of the stabilisation scale is sufficient to suppress higher dimensional condensates. Ioffe's second 
criticism is based on a calculation[25] where radiative corrections are not properly implemented and a 
number of other errors are made. The correct results are given in ref. [13], where we extend our previous 
analysis of the 'proton spin' problem systematically beyond the chiral limit using a new set of generalised 
Goldberger-Treiman relations. 
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one of the main advantages of the CPV method is that some non-perturbative information 
which is generic to QCD, i.e. independent of the target, is factored off into the composite 
operator propagator. This allows us to distinguish between non-perturbative mechanisms 
which are generic to all QCD processes and those which are specific to a particular target. 
As explained above, our contention is that the anomalous suppression in the first moment of 
g\ is of the first, target-independent, type. This conjecture could in principle be tested by 
DIS with non-nucleon targets, which may effectively be realised in semi-inclusive polarised 
DIS. This will be discussed in section 6. 

Both the parton and CPV methods allow a natural conjecture in which the origin of 
this suppression is attributed to 'glue' - either through a large polarised gluon distribution 
Ag(Q 2 ) in the parton description or due to an anomalous suppression of the slope of 
the topological susceptibility x'(0) in the CPV description. These conjectures are based 
on assumptions that the appropriate RG invariant quantities, Aq s or r$ 5PP , obey the 
OZI rule. The motivation for this is particularly strong in the CPV case, where it is 
supported by a range of evidence from low-energy meson phenomenology in the Ua(X) 
channel. Moreover, it identifies a fundamental physical mechanism as being at the origin 
of the 'proton spin' suppression - the screening of topological charge by the QCD vacuum. 

The two approaches therefore provide related, but complementary, insights into the 
nature of the 'proton spin' effect. Clearly, both insights are needed and both methods have 
a full part to play in understanding this intriguing and subtle phenomenon. 



5. Experiment 

In the last year, the SMC collaboration have completed their analysis of the final data 
from the 1996 run. This is shown in Fig. 6. 

gf at measured Q 2 preliminary 

O 93 
□ 96 

A 93+96 combined 



10 2 10 1 X 1 

Fig.6 SMC data for g{ 

Including an uncontroversial extrapolation in the unmeasured region 0.7 < x < 1, 
and taking into account the RG scaling from the measured Q 2 for a particular x to the 
reference value Q 2 = 10 GeV 2 , SMC quote[28] the following result for the first moment of 



4 

i 

2 
1 
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g\ in the measured range x > 0.003: 

T\(Q 2 = 10GeV 2 )l 00Q3 = / dx gf(x;Q 2 ) = 0.141 ±0.006 ±0.008 ±0.006 (5.1) 

io.003 

where the first error is statistical, the second is systematic and the third is due to the 
uncertainty in the Q 2 evolution. 




-5 ~ -4 -3 " -2 -1 

10 10 10 10 10 X 



Fig.7 SMC data for xg± including 'Regge' and 'QCD' small X extrapolations. 

The result for the entire first moment depends on how the extrapolation to the un- 
measured small x region x < 0.003 is performed. See Fig. 7. This is still a controversial 
issue. Using a simple Regge fit, SMC find T\ = 0.142 ± 0.017 from which they deduce 
a = 0.34 ± 0.17. Alternatively, using a small x fit in which the parton densities are 
parametrised at small x at a low Q 2 scale and then extrapolated to higher Q 2 using 
the perturbative QCD evolution equations [16, 29], SMC quote = 0.130 ± 0.017 and 
a = 0.22 ±0.17 (all at Q 2 = lOGeV 2 ). The key feature of the evolution equations is that 
while AS(x, Q 2 ) falls with increasing Q 2 , the polarised gluon distribution Ag(x, Q 2 ) rises. 
The net effect is that g\{x, Q 2 ) is driven strongly negative at Q 2 = lOGeV 2 for sufficiently 
small x. This gives a potentially large negative contribution to the first moment T^, but 
with relatively large errors[29]. 

Moreover, SMC have recently published an alternative analysis [30] of their data, using 
a different event selection, this time quoting a slightly lower number for the integral over 
the measured region of x: 

T\(Q 2 = 10GeV 2 )l 00Q3 = f dx g{{x;Q 2 ) = 0.133 ±0.005 ±0.006 ±0.004 (5.2) 

io.003 

Clearly it is premature to draw too strong a conclusion given the large errors on the 
experimental determinations of and a and the uncertainty over the small x extrapola- 
tion. Future studies of the small x region of g\ and the Q 2 scaling behaviour of the gluon 
distribution Ag(x, Q 2 ) will be important challenges to experimentalists. Nevertheless, it is 
extremely encouraging that the CPV prediction (4.15) is firmly in the region favoured by 
the data. This gives us extra confidence that the explanation of the 'proton spin' problem 
in terms of topological charge screening by the QCD vacuum is correct. 
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6. Semi-Inclusive Polarised DIS 



Recently, a new proposal to exploit semi-inclusive DIS in the target fragmentation 
region to elucidate the 'proton spin' effect has been presented [10, 11]. The idea is to test 
the mechanism of topological charge screening, or more precisely the prediction of 'target 
independence', suggested by the CPV method by using semi-inclusive DIS in effect to make 
measurements of the polarised structure functions of other hadronic targets besides the 
proton and neutron. 

The essence of the target independence conjecture is that for any hadron, the singlet 
axial charge in eq.(2.3) can be substituted by its OZI value multiplied by a universal 
(target-independent) suppression factor s(Q 2 ) determined, up to radiative corrections, by 
the anomalous suppression of the first moment of the topological susceptibility a/ x'(0). 
For example, for a hadron containing only u and d quarks, the OZI relation is simply 
a = a 8 , so we predict: 

ri = ^r(a 3 + ^(l + 4,)a 8 ) (6.1) 

where 

S[Q ] ~ C? s (a s ) a 8 {b ' 2) 

Since s is target independent, we can use the value measured for the proton to deduce Ti 
for any other hadron simply from the flavour non-singlet axial charges. From our spectral 
sum rule estimate of a/x'(0), we find s ~ 0.6 at Q 2 = 10 GeV 2 , while the central values of 
the results for a°(Q 2 ) extracted by SMC from the data give even lower values for s. 

The non-singlet axial charges for a hadron B are given by the matrix elements of 
the flavour octet axial currents, so can be factorised into products of SU(3) Clebsch- 
Gordon coefficients and reduced matrix elements. Together with the target independence 
conjecture, this allows predictions to be made for ratios of the first moments of the polarised 
structure functions gf for different £>, which involve only group theoretic numbers and the 
universal suppression factor s. Some of the most interesting are: 

= 2.-l-3(2. + l)F/Z> 
11 1 2s + 2-QsF/D V ; 

rr/rf = r?>? = |±| (6.4) 

where (E°) is the state with valence quarks uuc (ddc). The results for A and E c are 

particularly striking because of the 2s — 1 denominator factor, which is very small for the 
range of s favoured by experiment. These examples therefore show spectacular deviations 
from the valence quark counting (OZI) expectations, which would give the ratio 4. They 
also turn out to be the ones which allow the most clear-cut experimental interpretation. 

The proposal of ref.[10] is that these ratios can be realised in semi-inclusive DIS in a 
kinematical region where the detected hadron h (a pion or D meson in these examples) 
carries a large target energy fraction, i.e. z approaching 1, with a small invariant momentum 
transfer. To understand this, we briefly review some of the theory of semi-inclusive DIS. 
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Fig. 8 Semi-inclusive DIS: current fragmentation region. 




Fig. 9 Semi-inclusive DIS: target fragmentation region. 

There are two distinct contributions to the semi-inclusive DIS reaction eN — > ehX, 
coming from the current and target fragmentation regions. These are shown in Figs. 8 
and 9. The current fragmentation events are described by parton fragmentation func- 
tions D^(y3^, Q 2 ), where i denotes the parton, while the target fragmentation events are 
described by fracture functions[31] Mj lN (x, z, Q 2 ) representing the joint probability distri- 
bution for producing a parton with momentum fraction x and a detected hadron h (with 
momentum p' 2 ) carrying energy fraction z from a nucleon N. The lowest order cross section 
for polarised semi-inclusive DIS is: 



x 



dAa Y P 47ra 2 



dxdydz 



i 



AM^ N {x,z, Q 2 ) + 



1-x 



Aq t {x,Q 2 )D 



1 — x 



(6.5) 



where z = Eh/ En (in the photon-nucleon CM frame). The notation is slightly different 
from section 3. Here, Aqi(x) refers to quarks and antiquarks separately and a sum over 
both is implied. AM^fi, 5, Q 2 ) is the polarisation asymmetry of the fracture function. 
The NLO corrections to eq.(6.5) are given in ref.[32]. 

In fact, for our purposes it is better to use the extended fracture functions introduced 
recently in refs.[33]. These new fracture functions AM- bN (x, z, t, Q 2 ) have an explicit 
dependence on the invariant momentum transfer t. The original fracture functions are 
found by integrating over t in the range t < 0(Q 2 ). One of the advantages of the extended 
fracture functions is that they obey a simple, homogeneous, RG evolution equation: 



duj 



APij , a s ) AMj N (u, x, t, Q 2 ) (6.6) 
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The differential cross section in the target fragmentation can then be written, analogously 
to eq.(2.1), as 

dAa target y 4 2 

X ~dxdydzdt ~ ~2 s ^ ( x ' 2 ' ' *3 ) 

where AM t is the fracture function equivalent of the inclusive structure function , 
(c.f. eq.(3.2)): 

AMr( V) t,Q 2 ) = ^e 2 AM^(x,M,Q 2 ) (6.8) 

i 

For ^ approaching 1, i.e. the hadron carrying a large target energy fraction, this target 
fragmentation process may be simply modelled by a single Reggeon exchange (see Fig. 10). 
This corresponds to the approximation 

AM? N (x,z,t,Q 2 ) ~ F(t)(l-z)- 2 ^gf(x B ^Q 2 ) (6.9) 

The notation here is z = p'i-llVi-1-, %B = Q 2 /2k.q, l — z = x/x&, and i = — /c 2 <C 0(Q 2 ) so 
that ^ ~ 5. The Reggeon emission factor F(t)(l — z)~ 2oiB ^ cancels in the ratios of cross 
sections we consider. 



e(Pi) 




P (P 2 ) 

h(pj) 

Fig. 10 Single Reggeon exchange model of ep — > ehX . 

Although single Reggeon exchange is only an approximation to the more fundamental 
QCD description in terms of fracture functions, it shows particularly clearly how observing 
semi-inclusive processes at large z with particular choices of h and iV amounts in effect 
to performing inclusive DIS on virtual hadronic targets B. Since the predictions (6.3,6.4) 
depend only on the SU(3) properties of £>, together with target independence, they will 
hold equally well when B is interpreted as a Reggeon rather than a pure hadron state. For 
example, from the quark diagram in Fig. 11, we easily see that the semi-inclusive reaction 
ep — > en~ X measures gf for B — A ++ . 

u )a++(10) 

u J 




Fig. 11 Quark diagram for the NhB vertex in the reaction ep — > en X where B has the quantum 
numbers of A++. 
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We can therefore make simple predictions for the ratios 7Z of the first moments of the 
polarised fracture functions J Q 2 dxAM± N (x, z, t, Q 2 ) ~ Jq dxBgf(xB,t,Q 2 ) for various 
reactions. We emphasise again that these do not depend on any detailed model of the 
fracture functions, such as single or even multi-Reggeon exchange. (A more formal justi- 
fication in terms of the fracture function - cut vertex equivalence [33] is currently under 
study.) The ratios (6.3,6.4) are obtained in the limit as z approaches 1, where the reaction 
eN — > ehX is dominated by the process in which most of the target energy is carried 
through into the final state h by a single quark. We therefore predict 



ep — > en X\ 



2s + 2 



en 



en+X J z~i 2s - 1 



(6.10) 



Similarly for charmed mesons: 



K 



ep 



en 



eD~X \ 
■ eD°X ) 



2s + 2 
z~i 2s - 1 



(6.11) 



For strange mesons, on the other hand, the ratio depends on whether the exchanged 
object has SU(3) quantum numbers in the 8 or 10 representation, so the prediction is less 
conclusive, e.g. 



K 



ep 



zK°X\ 



en 



eK+X J z~i 



2s-l-3(2s + l)F*/D* 



2s -1 
2s + 1 
2s -1 



3(2s 
(10) 



1)F*/D* 



(8) 



(6.12) 



At the opposite extreme, for z approaching 0, the detected hadron carries only a small 
fraction of the target nucleon energy. In this limit, the ratio of cross section moments for 
ep — > en~(D~)X and en — > en + (D°)X is simply the ratio of the structure function 
moments for the proton and neutron. 

Interpolating between these limits, we expect the ratio 1Z in the range < z < 1 for 
the reactions en — > en + (D°)X over ep — > eTc~(D~)X to look like the sketch in Fig. 12. 



0.30 
0.20 
0.10 





-0.10 
-0.13 
-0.20 
-0.28 
-0.30 



Fig. 12 Cross section ratios 1Z for en — > e7T + (D°)X over ep 
Z — > 1, contrasting the OZI and CPV predictions 




en {D )X between Z — > and 
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The difference between the OZI (or valence quark model) expectations and these 
predictions based on our target-independent interpretation of the 'proton spin' data is 
therefore quite dramatic, and should give a clear experimental signal. 

Since the proposed experiments require particle identification in the target fragmen- 
tation region, they are difficult to do at a polarised fixed-target experiments such as COM- 
PASS, or even HERMES, which are better suited to studying semi-inclusive processes in 
the current fragmentation region. A better option would be a polarised ep collider, such 
as HERA[34]. To conclude, we shall make some brief comments about the experimental 
requirements necessary for these predictions to be tested[ll]. 

The first requirement is to measure particles at extremely small angles (9 < 1 mrad) , 
corresponding to t less than around 1 GeV 2 . This has already been achieved at HERA 
in measurements of diffractive and leading proton/neutron scattering. The technique for 
measuring charged particles involves placing detectors commonly known (but see the dis- 
cussion session) as 'Roman Pots' inside the beam pipe itself. In the present experiments, 
this detection system is known as the Leading Proton Spectrometer (LPS). 

The next point is to notice that the considerations above apply equally to p as to it 
production, since the ratios 1Z are determined by flavour quantum numbers alone. The 
particle identification (ID) requirements will therefore be less stringent, especially as the 
production of leading strange mesons from protons or neutrons is strongly suppressed. 
However, we require the forward detectors to have good acceptance for both positive and 
negatively charged mesons M = n, p in order to measure the ratio (6.10). 

The reactions with a neutron target can be measured if the polarised proton beam 
is replaced by polarised 3 He. In this case, if we assume that 3 He = Ap + Bn, the cross 
section for the production of positive hadrons h + measured in the LPS is given by 

a( 3 He -> h + ) ~Aa(p^ h + ) + Ba{n -> p) + Ba(n -> M+) (6.13) 

The first contribution can be obtained from measurements with the proton beam. How- 
ever, to subtract the second one, the detectors must have sufficient particle ID at least to 
distinguish protons from positively charged mesons. 

Finally, estimates of the total rates suggest that around 1% of the total DIS events 
will contain a leading meson in the target fragmentation region where the LPS has non- 
vanishing acceptance (z > 0.6) and in the dominant domain x < 0.1. The relevant cross 
sections are therefore sufficient to allow the ratios 1Z to be measured. 

We conclude that while these proposals undoubtedly pose a challenge to experimen- 
talists, they are nevertheless possible. Given the theoretical importance of the 'proton 
spin' problem and the topological charge screening mechanism, there is therefore strong 
motivation to perform target fragmentation experiments at polarised HERA. 
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